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Hybrid functionals for solids with an optimized Hartree-Fock mixing parameter 
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(Screened) hybrid functionals are being used more and more for solid-state calculations. Usually 
the fraction a of Hartree-Fock exchange is kept fixed during the calculation, however there is no 
single (universal) value for a which systematically leads to satisfying accuracy. Instead, one could 
use a property of the system under consideration to determine a and in this way the functional 
would be more flexible and potentially more accurate. Recently, it was proposed to use the static 
dielectric constant e for the calculation of a [Shimazaki and Asai, Chem. Phys. Lett. 466, 91 (2008) 
and Marques et al., Phys. Rev. B 83, 035119 (2011)]. We explore this idea further and propose 
a scheme where the connection between e and a is optimized based on experimental band gaps. 
e, and thus a, is recalculated at each iteration of the self-consistent procedure. We present results 
for the band gap and lattice constant of various semiconductors and insulators with this procedure. 
In addition, we show that this approach can also be combined with a non-self-consistent hybrid 
approximation to speed up the calculations considerably, while retaining an excellent accuracy in 
most cases. 

PACS numbers: 71. 15. Mb, 71.20.-b, 77.22.Ch 
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I. INTRODUCTION 

Hybrid functionals^, and in particular the screened 
versions^— provide, in combination with a generaliza- 
tion (Ref. 6) of the Kohn-Sham equations^ an efficient 
way to deal with the well-known band gap problem^ in 
density functional theory (DFT) J£ In hybrid functionals, 
(semi)local [i.e., local density (LDA) or generalized gra- 
dient approximation (GGA)] and Hartree-Fock (HF) ex- 
change are mixed, and since the trends of pure semilocal 
and pure HF methods is to underestimate and overesti- 
mate band gaps, respectively, the mixing of both usually 
leads to more accurate band gaps. The use of hybrid 
functionals has been justified by (more or less) formal 
arguments (Refs. Q] and 111] ) . 

One well-established screened hybrid functional for 
solids is the one proposed by Heyd, Scuseria, and 
Ernzerhof 4 (HSE), which is based on the GGA functional 
of Perdew, Burke, and Ernzerhof 12 (PBE): 



£U = aE. 



SR,HF 



■(l-a)E: 



SR.PBE 
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LR.PBE 



E, 



PBE 



(1) 



where a (g [0, 1]) is the mixing parameter for the short- 
range exchange (in screened hybrid functionals the long- 
range exchange is pure semilocal). The indices "x" and 
"c" denote the exchange and correlation energy contribu- 
tions. The HSE functional is usually used with a = 0.25. 
Another parameter in Eq. {T]) is the screening parameter 
A in the error function (A = 0.11 bohr -1 in the ver- 
sion HSEO&i^), which determines the separation between 
short-range and long-range exchange. A larger amount 
a of HF exchange increases the band gap, while a larger 
value for A (i.e., more screening) decreases it. Actually, 
this means that depending on the choice of these two 



parameters more or less any desired result between PBE 
and HF can be obtained. Furthermore, it also means that 
a reasonable way to fix the parameters is very important. 
In the HSE06 functional, the first of these parameters (a) 
was chosen on the basis of theoretical considerations?^ 
whereas the second one (A) was fitted to experimental 
results^ It has been shown that by fitting both parame- 
ters further improvement can be achieved^ 

At this point, it is worth recalling that already in 1990 
Bylander and Klcinman 3 proposed an LDA-based func- 
tional containing screened HF exchange (sX-LDA): 
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which can be regarded as a hybrid functional with a — 
1. In the sX-LDA functional, the long-range and short- 
range exchange are split using the exponential function 
(Yukawa potential) . 

Up to now, the majority of calculations with hybrid 
functionals (full-range, short-range, or long-range) have 
been done with fixed values for a and A. However, al- 
lowing a or A to depend on a property of the system is a 
way to make the functional more flexible and thus poten- 
tially more accurate. A brief summary of such schemes is 
mentioned below. The sX-LDA functional is usually used 
with a screening parameter A which is calculated using 
the average of the valence electron density (see, e.g., Ref. 
\la for recent calculations). In Ref. [lg, it was proposed 
to make a position-dependent by using the electron den- 
sity p, its derivative Vp, and the kinetic-energy density, 
while in Ref. Il7l the use of a position-dependent A which 
depends on p and Vp has been proposed. Shimazaki 
and Asaii£~— proposed several functionals in which ei- 
ther both a and A or only a are determined using the 



static dielectric constant e s . In the method presented in 
Ref. 21 by Stein et al, the screening parameter is tuned 
such that the Koopmans' theorem (which requires the 
use of orbital eigenvalues and total energies) is obeyed 
as closely as possible. Marques et al£2- considered two 
ways for the calculation of a: either with the static di- 
electric constant e s or with the average of |V/o| j p in the 
unit cell2£ (as done originally in a similar context for the 
modified Becke- Johnson exchange potential 2 ^) . We also 
mention that interesting discussions about the link be- 
tween screened hybrid functionals and quasiparticle the- 
ories can be found in Refs. [22J andH. 

In the present work, we further explore the use of 
the dielectric constant for the calculation of the fraction 
of HF exchange in the screened hybrid functional YS- 
PBEO,— which is based on the Yukawa operator. We will 
show the results obtained with this scheme for the band 
gap and lattice constant of solids. Our work is organized 
as follows. In Sec. |TT]the computational details are given, 
while the description of our method and the results are 
presented in Sec. IIIII Finally, the summary of our work 
will be given in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

All calculations were performed with the WIEN2k 
software,— which is based on the full potential (lin- 
earized) augmented plane wave and local orbitals 
method 25 for quantum calculations on periodic systems. 
Recently, unscreened and screened hybrid functionals 
were implemented into WIEN2k4 In screened hybrid 
functionals, the screening of the Coulomb operator is 
done by using the exponential function, and it was shown 
(Ref. M) that by choosing carefully the screening param- 
eter A, the results are very close to the results from the 
error function-based screened hybrid functionals (e.g., 
HSE06). More specifically, the screening parameter used 
in the exponential function should be about 1.5 times 
larger than the one used for the error function (see Refs. 
I5l ll8l - l20l for details) . The results presented in the present 
work were obtained with the YS-PBEO functional which 
is based on the PBE functional^ 2 - and was used in our pre- 
vious worksi 26 i 27 The YS-PBEO calculations were done 
with a fixed screening parameter of A = 0.165 bohr -1 , 
which gives results close to the results obtained with the 
HSE06 functional^ Spin-orbit coupling has not been 
considered explicitly in hybrid functional calculations, 
but its influence on the band gap is considered to be 
the same as for the PBE functional and the gaps have 
been corrected accordingly. 

The imaginary part of the dielectric function e was cal- 
culated using Fermi's golden rule and the independent 
particle approximation. The Kramers-Kronig transfor- 
mation was used to obtain the real part of e. Details of 
this approach can be found in Ref. [28|. This gives the 
G = G' = element of the dielectric matrix from which 
the static macroscopic dielectric constant (which is a ten- 



sor since solids are anisotropic) can be obtained by taking 
the limit uj — > and q — > of the real part. For solids 
whose symmetry leads to vanishing off-diagonal elements 
of this tensor we take the geometric mean of the diagonal 
elements of this tensor and from now on we will refer to it 
simply as the dielectric constant and use the symbol e* . 
Systems with non-vanishing off-diagonal elements were 
not considered in the present work. All parameters of 
the calculations, such as basis-set size or Brillouin zone 
sampling were tested for convergence. Except for the cal- 
culations of the lattice constants, the experimental struc- 
tures have been used. 



III. RESULTS AND DISCUSSION 



A. Band gap 
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FIG. 1: (Color online) The fraction of HF exchange required 
to reproduce the experimental band gap (a exp ) versus the 
inverse of the dielectric constant (1/e*) which is obtained from 
a YS-PBEO calculation using the corresponding a cxp . The 
straight line [a op t = o op t(e*)] is a fit corresponding to the 
least square error of the gap and not of a. 

Hybrid functionals with fixed amount of exact ex- 
change and screening length can be problematic if ap- 
plied to different types of systems. For instance HSE 
is knownsi to work well for many small-gap semicon- 
ductors but strongly underestimates the gaps of highly 
ionic compounds such as NaCl. This means that in or- 
der to find an hybrid functional with broader applica- 
bility, it is important to include different types of solids 
in the fitting set. The solids that we considered for this 
study are listed in TablelU They include highly ionic com- 



TABLE I: Fundamental band gaps (in eV) of 24 solids categorized into four different groups: ionic compounds, sp- 
semiconductors, transition- metal oxides (labeled as "TmO"), and other transition- metal compounds (labeled as "TmX"). 
The Strukturbericht symbols are indicated in parenthesis. The three columns "YS-PBEO" show the results obtained with the 
hybrid functional YS-PBEO with a = 0.25 (Sec. IIII A|) . a opt (Sec. IHI A|) . and a opt with the non self-consistent diagonal-HF 
approximation (Sec. IIIIC[) . The experimental band gaps (see Ref. 29j for the references) are shown for comparison. In addi- 
tion, a op t (obtained from the fitting procedure), and the static dielectric functions obtained from YS-PBEO (a op t), PBE, and 
experiment are also shown. 
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"It is not recommended to apply the diagonal-HF approximation 
to ScN since it is metallic with PBE. 
^Reference 2g. 
^Reference l39l . 
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pounds with large band gaps, technologically important 
sp-semiconductors and more or less strongly correlated 
transition-metal compounds which can be split further 
into oxides (TmO) and non-oxides (TmX). 

As already mentioned, the dielectric constant e* will 
be used for the determination of the fraction a of HF 
exchange. Note that in the present work, the screening 
parameter A has been kept fixed at A = 0.165 bohr -1 
(Ref. |fj|) to make comparison with YS-PBEO (and there- 
fore HSE06) possible. As we can see from Fig. [TJ choos- 



ing £* seems to be a very good choice. Indeed, there 
is a nice linear correlation between ev. eX p, which is the 
amount of HF exchange required to reproduce the ex- 
perimental band gap and 1/e* which is obtained from a 
YS-PBEO calculation using the corresponding a cxp . We 
considered the linear relation a op t — A + B/e* and a fit 
procedure minimizing the least-square error in the gaps 
leads to A = 0.147 and B = 0.634. Note, that our lin- 
ear approximation does not go through the origin and 
the slope is not one, while in a previous work reported 
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FIG. 2: (Color online) Difference between calculated (full 
symbols: YS-PBEO with a op t, open symbols: YS-PBEO with 
a — 0.25) and experimental band gaps. The shape of the 
symbol indicates the type of solids. 



in Ref. 22 a strict proportionality a = l/e PBE was used 
(in unscreened PBEO). The linear fit is shown in Fig. [TJ 
where we can see that most data points are quite close to 
it. The exceptions are the ionic chlorides and the noble 
gases (not shown explicitly) whose a C xp are larger than 
the values a op t obtained from the fit and rutile and stron- 
tium titanate, where a oxp is a bit smaller than the fitted 
a opt . Thus, it can be expected that the gaps of the chlo- 
rides will be underestimated, whereas the gaps of rutile 
and strontium titanate will be overestimated. 

The resulting band gaps arc listed in Table [J and are 
compared with the values obtained with the standard 
fixed value a = 0.25. The band gaps are also plotted 
in Figs. [2] and [3] for a better comparison. It can be 
clearly seen that in nearly all cases, improved band gaps 
are obtained with cv opt compared to a = 0.25, and in 
some cases the improvement is quite impressive. In par- 
ticular, the band gaps of the highly ionic compounds, 
which are known 4 ^ to be underestimated by about 3 eV 
in HSE, are reproduced much more accurately with the 
e*-dependent a opt . As expected from Fig. [TJ the use of 
a op t leads to underestimations (which are much smaller 
than when HSE is used) for the chlorides and overestima- 
tions in the case of rutile and strontium titanate. For the 
other solids the performance can be considered as excel- 
lent. The deviation from experiment is within a range of 
±0.3 eV which is quite good, since experimental errors 
and temperature effects need to be considered as well. 



FIG. 3: (Color online) Same as Fig. [2] but only for solids 
with band gaps smaller than 7 eV. 



The performance is even slightly better than that of the 
modified Bccke- Johnson potential with the parameters 
suggested in Ref. [23- 

Since the static dielectric constant e* is a central quan- 
tity in the described approach, it is important to compare 
the calculated values to the experimental results. Ta- 
ble [J shows the values obtained from calculations (with 
PBE and YS-PBE0 Qopt ) and experiment. We can see 
that YS-PBEO with a opt underestimates the experimen- 
tal values by about 30%, contrary to PBE which slightly 
overestimates them. This is not surprising, since gaps are 
underestimated by PBE and lower gaps lead to higher 
dielectric constants. In our calculations, however, the di- 
electric function is calculated in the independent particle 
approximation and does not include electron-hole inter- 
actions, which would lead to an increase in the dielectric 
constants^ 1 - An improved calculation for e* would require 
a re-parameterization of a op t- 

As a last remark in this section, we also mention that 
an alternative way to determine a would be to use 
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which is the average of |Vp| / p in the unit cell and has 
been used successfully for the modified Becke- Johnson 
potential . 23 i 29 i 42 i 43 However, as Fig. [4] shows, the situa- 
tion with g is quite different since there is no clear correla- 
tion between the values of a eX p and g. Different classes of 
materials (sp-semiconductors, ionic chlorides or fluorides, 



o 



u,/ - 


■ 


ionic 




z 












Q) 


Z 




T 


sp 




o 


01 


0,6 - 


• 


TmO 
TmX 




■ ■ 

7; 


Tl 

-n ■ 

■ 

■ 










O 








LiCI ■ o 


i - 


0,5 - 








■ 


TJ 










-, N 












2 3 










N 


3 O 




0,4 - 






CO 


o 










• 


O 








> CO 
~° O 


T 


!= 

o 




0,3 - 




CO 


T 


A H 










• 


d 








Ge 


CO 

o 


A- 




0,2 - 




MoS 2 « 


z 






1 ' ' 


| i i i | 


1 1 1 1 


I I I | I I I 


| 1 1 1 | 1 1 1 | 



0,04 



1 1,2 1,4 1,6 1,8 2 2,2 2,4 

g [a.u.] 

FIG. 4: (Color online) The fraction of HF exchange required 
to reproduce the experimental band gap (a cxp ) versus g (see 
text for definition). 



transition metal compounds) would require drastically 
different parameterizations. The transition-metal oxides 
show only small variations in ~g but require quite different 
a opt , which is of course very problematic if any kind of 
relation a op t = aopt(<?) is desired. This is in quite strong 
contrast to the findings in Ref. [22J, where they found a 
good relation between a and g. However, their set of 
solids was much more restricted in terms of transition- 
metal or highly ionic systems, and furthermore, their cal- 
culations were obtained with the pseudopotential plane 
wave method, whereas the calculations presented in the 
present work were obtained from an all-electron method, 
this difference leading certainely to different values of g. 
Therefore, we considered only the parameterization in 
terms of the dielectric constant. 



B. Lattice parameters 

Another important test for functionals is the total en- 
ergy, in particular how accurate equilibrium structural 
parameters or atomization energies can be described. 
Note, that a op t will change as function of the lattice pa- 
rameter and this could have important consequences on 
the equilibrium geometry. Therefore, we calculated the 
equilibrium lattice parameters for a few selected solids by 
calculating the total energies for different lattice param- 
eters and fit them to the Birch-Murnaghan 4 ^ equation 
of state (EOS). The minimum of this fit corresponds to 
the equilibrium lattice parameter. An example is shown 
in Fig. [5] for silicon and smooth total energy fits can 
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FIG. 5: (Color online) Total energy of silicon as a function 
of atomic volume for 3 different functionals. The lines rep- 
resent fits to the Birch-Murnaghan-EOS, the arrows indicate 
the corresponding minima. 



be obtained for all functionals, including YS-PBE0 with 
a op t- 

Table [TT] shows the equilibrium lattice parameters ob- 
tained with the hybrid functional YS-PBEO for a few 
representative solids. For comparison, the results ob- 
tained with other functionals, namely, LDA, Wu-Cohen 
(WC), 48 and HSE06, as well as from experiment are also 
shown. The selected solids include cases with low (e.g., 
LiCI), intermediate (e.g., ScN), and high (e.g., Si) values 
of the dielectric constant. As already knownj 4 ^— LDA 
underestimates strongly the lattice constants, while WC 
belongs to the group of the most accurate GGA func- 
tionals for this property, and HSE06 has the tendency 
to give too large values (albeit not as much as PBE). 
In most cases YS-PBEO with fixed a = 0.25 is in ex- 
cellent agreement with HSE06, while for NaCl a slightly 
larger deviation is observed. The results for YS-PBEO 
with a op t = A + B/e* are in most cases very similar to 
those with fixed a and thus one could conclude that a 
hybrid functional with optimized a is well suited also for 
lattice parameter determinations. However, we can see 
that for the highly ionic compounds LiCI and NaCl, a 
dramatic increase of the equilibrium lattice parameters, 
reaching completely unphysical values (an increase from 
5.15 to 5.55 A and 5.69 to 6.20 A for LiCI and NaCl, re- 
spectively), is obtained. Clearly this approach is not rec- 
ommendable in such a case. In order to know where the 
problem comes from it is necessary to look more closely at 



TABLE II: Comparison between theoretical and experimental lattice parameters (in A) for several solids. All experimental 
values except ScN are corrected for zero-point anharmonic expansion.— The values of a opt and da op t/d\n V at the experimental 
volume are also shown. 



Solid 



exp. 



LDA 



WC 



HSE06 



YS-PBEO 


YS-PBEO 


(a = 0.25) 


(a op t) 


3.55 


3.55 


5.46 


5.44 


4.36 


4.34 


3.61 


3.60 


4.52 


4.50 


4.51 


4.52 


4.82 


4.86 


5.15 


5.55 


5.69 


6.30 



*opt(V = Kxp) 



-{V = Vcxp) 



c 

Si 

SiC 

BN 

GaN 

ScN 

CaO 

LiCl 

NaCl 



3.54° 
5.42° 
4.34° 
3.59° 
4.52° 
4.50 c 
4.79" 
5.07" 
5.60" 



"Reference 
'Reference 
c Reference 



3.54° 
5.41° 
4.33" 
3.59" 
4.46° 
4.43 c 
4.72" 
4.97° 
5.48" 



3.56° 
5.44° 
4.36" 
3.61" 
4.50° 
4.47 c 
4.78° 
5.07" 
5.62° 



3.55" 
5.44* 
4.35* 
3.60 6 
4.49* 



5.12" 
5.64* 



0.306 
0.224 
0.284 
0.346 
0.323 
0.254 
0.427 
0.447 
0.498 



-0.025 

-0.015 

-0.041 

-0.021 

-0.098 

0.042 

0.077 

0.15 

0.13 



the volume dependence of a opt (see Table ILl]) . In the case 
of CaO and especially LiCl and NaCl there is a very large 
slope of a op t as function of the volume, indicating that 
e* (and the band gap) increases (decreases) strongly with 
reduced volume. On the other hand, for typical semicon- 
ductors cv op t varies much less strongly when the volume is 
changed and the variations in the band gap are compen- 
sated by those in the dielectric function. Together with 
the fact that a pt is almost twice as large as «hse = 0.25 
the equilibrium lattice parameters become very large in 
highly ionic materials. It should be noted, however, that 
it is the variation of a pt with volume and not the large 
a value itself, which causes the failure to obtain good 
equilibrium volumes. For instance, for LiCl a lattice pa- 
rameter of 5.17 A is obtained if the fixed value a = 0.447 
(obtained at V eX p) is used. Therefore, an approach using 
a fixed a offers an alternative way to describe simultane- 
ously the energy band gap and structural parameters in 
a reasonable way. 

Next we want to test the YS-PBEO (a op t) scheme for 
a more complicated example: the pressure-induced Bl- 
B2 phase transition in CaO. Here the goal is to find the 
transition pressure above which the enthalpy of the B2 
phase becomes lower than that of the Bl phase. Experi- 
mentally this transition is found to occur in the range of 
60-70 GPa4^ Standard DFT functionals can reproduce 
this quite well: LDA predicts the transition to occur at 
57 GPa, WC at 61 GPa and PBE is in the center of 
the experimental range with a transition pressure of 66 
GPa. Traditional YS-PBEO with a = 0.25 also predicts 
it correctly at 67 GPa, but if ev op t is used the transition 
pressure is shifted to 87 GPa. This means that the re- 
sult is worsened in a similar way as the result for the 
CaO lattice parameter (it is however not as bad as the 
0.5 A error for LiCl). The shift of the transition point to 
higher pressure is related to the different a op t values of 
both phases (Bl: 0.43, B2: 0.39). The higher a opt of the 



Bl-phase reduces its total energy so that more pressure 
is required for a transition to the B2 phase. 



C. Combining the optimized-a-approach with the 
diagonal- hybrid-approximation 

Since hybrid calculations are computationally demand- 
ing, it is desirable to have a way to get similar results 
to the approach described above in a much shorter time. 
This is made possible by the non self-consistent diagonal- 
only hybrid approximation proposed in Ref. l5Ct In 
this approximation, only the diagonal elements of the 
perturbation Hamiltonian (hybrid minus semilocal) are 
calculated (in the basis of the semilocal orbitals), while 
the non-diagonal elements are neglected. This saves the 
time of evaluating most of the matrix elements and also 
a self-consistency cycle is not necessary. It was found 
previously^ that for common semiconductors and insu- 
lators this procedure leads to gaps in very close agree- 
ment with fully self-consistent calculations. In addition 
we assume a linear dependency of the gap and 1/e* on a, 
so that we can obtain a opt and the corresponding band 
gap from two calculations at a\ — (plain PBE) and 
a 2 = 0.3. 

The results obtained by this procedure are included in 
Table HI Except for the transition metal oxides (and to 
a much lesser degree the ionic fluorides), the deviation 
from the full-hybrid gaps is quite small, usually in the 
range of a few hundredth of an eV. On the other hand, 
the reduction of the computational effort is substantial, 
which allows to apply this approximation also to large 
unit cells. 

In order to analyze why this scheme works for most 
cases but has some exceptions, we take a closer look at 
CdS as an example for a working case and at ZnO as an 
example for an exception. The relevant data are shown 
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difference between the dielectric constant obtained with 
full hybrid calculations and that obtained with the diag- 
onal approximation. The reason for that is that in the 
diagonal calculations only the eigenvalues change (to al- 
most the same values as in full hybrid calculations), but 
the orbitals do not and still belong to the single-particle 
Hamiltonian with a semilocal potential. While in the case 
of CdS the eigenfunctions from PBE or full hybrid-DFT 
calculations do not differ much so that the dielectric con- 
stants with these two schemes remain very similar, in the 
case of ZnO there is a strong difference in the eigenfunc- 
tions from PBE or hybrid-DFT. This leads to different 
momentum matrix elements and therefore dielectric con- 
stants and finally the band gaps differ substantially. 



IV. SUMMARY AND CONCLUSION 



FIG. 6: (Color online) Dependence of the band gap and the 
inverse dielectric constant on the amount of HF exchange in 
full and diagonal YS-PBEO in CdS. Lines are a guide to the 
eye. 
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FIG. 7: (Color online) Dependence of the band gap and the 
inverse dielectric constant on the amount of HF exchange in 
full and diagonal YS-PBEO in ZnO. Lines are a guide to the 
eye. 



A hybrid functional was presented which, in contrast 
to most currently established hybrid functionals, uses a 
HF mixing parameter which is individually adapted for 
each investigated system. This is achieved automatically 
using the calculated dielectric constant. Since this is a 
global quantity of a bulk material, the functional as de- 
scribed here can be applied reasonably only to systems 
with strict 3-dimensional periodicity. However it may 
be generalized in the course of future work so that de- 
fect structures, molecules, surfaces and interfaces can be 
dealt with. Actually, it has been claimed^ that a hybrid 
functional with position-independent mixing factor and 
screening length is not appropriate for such cases. 

The results obtained with this functional are quite im- 
pressive. Band gaps for a wide range of materials are 
much better reproduced than by hybrid functionals with 
a fixed parameter. Structural properties, which are usu- 
ally not considered in similar studies, are in most cases of 
excellent quality and only in highly ionic insulators prob- 
lems appear. However, even in these cases, good results 
can be obtained when one fixes a op t at one volume and 
neglects its volume variations. 

We have also tested an approximate version of this 
approach, namely a non self-consistent hybrid scheme 
considering only the diagonal HF matrix elements. Ex- 
cept for transition metal oxides (where even non self- 
consistent GW methods fail) this approach leads to al- 
most identical band gaps and such a method can be used 
to predict reliable band gaps in systems with hundreds 
of atoms. 



in Figs. [6] and [7J We can see that in both cases the de- 
pendency of the band gaps on a using the full hybrid or 
the diagonal-hybrid approximation is nearly identical and 
fairly linear. For CdS the same holds for the dependency 
of 1/e* on a. However, for ZnO there is a substantial 
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